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D ■ Abstract 

The differential cross-section for the inclusive production of ip{2S) mesons in pp 
ps| I collisions at y/s=7TeV has been measured with the LHCb detector. The data 

• sample corresponds to an integrated luminosity of 36pb~^. The ip{2S) mesons are 

reconstructed in the decay channels ip{2S) — t- z^"^^" and ip{2S) — t- J/^/'tt^vt", with 
fSJ I the J/ip meson decaying into two muons. Results are presented both for promptly 

produced ip{2S) mesons and for those originating from 6-hadron decays. In the 
^ . kinematic range pt{iP{2S)) < 16GeV/c and 2 < y{ip{2S)) < 4.5 we measure 

^' aprompt(V'(25)) = 1.44 ± 0.01 (stat)± 0.12 (syst)l°iO (pol) ^b, 

abi'ip{2S)) = 0.25 ± 0.01 (stat) ± 0.02 (syst) ^b, 

^ ' where the last uncertainty on the prompt cross-section is due to the unknown ^(25) 

, polarization. Recent QCD calculations are found to be in good agreement with our 

measurements. Combining the present result with the LHCb J/ip measurements we 
determine the inclusive branching fraction 



B{b ij{2S)X) = (2.73 ± 0.06 (stat) ± 0.16 (syst) ± 0.24 (BF)) x 10" 
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where the last uncertainty is due to the B{b — )■ J/ipX), 13{J/ip — )■ ) and 

B{il){2S) — )• e^e~) branching fraction uncertainties. 
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1 Introduction 



Since its discovery, heavy quarkonium has been one of the most important test laboratories 
for the development of QCD at the border between the perturbative and non-perturbative 
regimes, resulting in the formulation of the nonrelativistic QCD (NRQCD) factorisation 
formalism [HEj. However, prompt production studies carried out at the Tevatron collider 
in the early 1990s [3] made clear that NRQCD calculations, based on the leading-order 
(LO) colour-singlet model (CSM), failed to describe the absolute value and the transverse 
momentum (px) dependence of the charmonium production cross-section and polarization 
data. Subsequently, the inclusion of colour-octet amplitudes in the NRQCD model has 
reduced the discrepancy between theory and experiment, albeit at the price of tuning 
ad hoc some matrix elements |2]. On the other hand, recent computations of the next-to- 
leading-order (NLO) and next-to-next-to-leading-order (NNLO) terms in the CSM yielded 
predictions in better agreement with experimental data, thus resurrecting interest in the 
colour-singlet framework. Other models have been proposed and it is important to test 
them in the LHC energy regime [HIS]. 

Heavy quarkonium is also produced from 6-hadron decays. It can be distinguished 
from promptly produced quarkonium exploiting its finite decay time. QCD predictions 
are based on the Fixed-Order-Next-to-Leading-Log (FONLL) approximation for the hh 
production cross-section. The FONLL approach improves NLO results by resumming px 
logarithms up to the next-to-leading order [HIEI. 

To allow a comparison with theory, promptly produced quarkonia should be separated 
from those coming from 6-hadron decays and from those cascading from higher mass states 
(feed-down). The latter contribution strongly affects J/ip production and complicates 
the interpretation of prompt J/ip data. On the other hand, ipi^S) charmonium has no 
appreciable feed-down from higher mass states and therefore the results can be directly 
compared with the theoretical predictions, making it an ideal laboratory for QCD studies. 

This paper presents a measurement of the ip{2S) meson production cross-section in 
pp collisions at the centre-of-mass energy ^/s = 7TeV. The data were collected by the 
LHCb experiment in 2010 and correspond to an integrated luminosity of 35.9±1.3 pb~^. 
The analysis is similar to that described in Ref. [S] for the J/tp production studies; in 
particular, the separation between promptly produced ip{2S) and those originating from 
6-hadron decays is based on the reconstructed decay vertex information. Two decay 
modes of the ip{2S) meson have been used: ^^{23) — )■ /i"*"//" and ^^{23) — t- J/ipn^n" 
followed by J/ip — >■ The J/ipii~^n~ mode, despite a larger background and a lower 

reconstruction efficiency, is used to cross-check and average the results, and to extend the 
accessible phase space. The production of ip(23) meson at the LHC has also been studied 
at the CMS experiment [9]. 

2 The LHCb detector and data sample 

The LHCb detector is a forward spectrometer [10], designed for precision studies of CP 
violation and rare decays of b- and c-hadrons. Its tracking acceptance covers approxi- 
mately the pseudorapidity region 2 < r/ < 5. The detector elements are placed along 
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the beam line of the LHC starting with the vertex detector, a sihcon strip device that 
surrounds the pp interaction region and is positioned at 8 mm from the beams during 
colhsions. It provides precise measurements of the positions of the primary pp interaction 
vertices and decay vertices of long-hved hadrons, and contributes to the measurement of 
particle momenta. Other detectors used for momentum measurement include a large area 
silicon strip detector located before a dipole magnet of approximately 4 Tm, and a com- 
bination of silicon strip detectors and straw drift chambers placed downstream. Two ring 
imaging Cherenkov detectors are used to identify charged hadrons. Further downstream 
an electromagnetic calorimeter is used for photon and electron detection, followed by a 
hadron calorimeter. The muon detection consists of five muon stations equipped with 
multi-wire proportional chambers, with the exception of the centre of the first station 
using triple-GEM detectors. 

The LHCb trigger system consists of a hardware level, based on information from 
the calorimeter and the muon systems and designed to reduce the frequency of accepted 
events to a maximum of 1 MHz, followed by a software level which applies a full event 
reconstruction. In the first stage of the software trigger a partial event reconstruction is 
performed. The second stage performs a full event reconstruction to further enhance the 
signal purity. 

The analysis uses events selected by single muon or dimuon triggers. The hardware 
trigger requires one muon candidate with a larger than 1.4GeV/c or two muon candi- 
dates with a pt larger than 560MeV/c and 480MeV/c. In the first stage of the software 
trigger, either of the two following selections is required. The first selection confirms the 
single muon trigger candidate and applies a harder cut on the muon at 1.8GeV/c. 
The second selection confirms the dimuon trigger candidate by requiring the opposite 
charge of the two muons and adds a requirement to the dimuon mass to be greater than 
2.5GeV/c^. In the second stage of the software trigger, two selections are used for the 
ip{2S) — )■ fi^fj^~ mode. The first tightens the requirement on the dimuon mass to be 
greater than 2.9GeV/c^ and it applies to the firtst 8pb~ of the data sample. Since this 
selection was subsequently prescaled by a factor five, for the largest fraction of the re- 
maining data (28pb~^) a different selection is used, which in addition requires a good 
quality primary vertex and tracks for the dimuon system. For the J/ip 7r~^7r~ mode only 
one selection is used which requires the combined dimuon mass to be in a ±120MeV/c^ 
mass window around the nominal J/tp mass. To avoid that a few events with high oc- 
cupancy dominate the software trigger CPU time, a set of global event cuts is applied 
on the hit multiplicity of each subdetector used by the pattern recognition algorithms, 
effectively rejecting events with a large number of pile-up interactions. 

The simulation samples used for this analysis are based on the Pythia 6.4 genera- 
tor [11] configured with the parameters detailed in Ref. [12]. The prompt charmonium 
production processes activated in Pythia are those from the leading-order colour-singlet 
and colour-octet mechanisms. Their implementation and the parameters used are de- 
scribed in detail in Ref. [12]. The EvtGen package [2] is used to generate hadron 
decays and the Geant4 package [13] for the detector simulation. The QED radiative 
corrections to the decays are generated using the Photos package [T6] . 
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3 Signal yield 



The two modes, ip{2S) — >■ fi~^^~ and ip(2S) — >■ J/t/^tt+tt", have different decay and back- 
ground characteristics, therefore dedicated selection criteria have been adopted. The 
optimisation of the cuts has been performed using the simulation. A common require- 
ment is that the tracks, reconstructed in the full tracking system and passing the trigger 
requirements, must be of good quality {x^/ndf < 4, where ndf is the number of degrees of 
freedom) and share the same vertex with fit probability -P(x^) > 0.5% {iIj{2S) — > 
and -P(x^) > 5% {ip{2S) — )■ J/ipn^Ti~). A cut > 1.2GeV/c is applied for the muons 
from the ip{2S) — )■ decay. For muons from J/tp{fi~^fi~)'n'~^'K~ we require a momen- 

tum larger than 8GeV/c and > 0.7GeV/c. Finally the rapidity of the reconstructed 
ipi^S) is required to satisfy the requirement 2 < y < 4.5. 

The ip(2S) —7- fi'^fi~ invariant mass spectrum for all selected candidates is shown 
in Fig. [U^a). The fitting function is a Crystal Ball [T7] describing the signal plus an 
exponential function for the background. In total 90600±690 signal candidates are found 
in the px range 0-12GeV/c. The mass resolution is 16.01±0.12 MeV/c^ and the Crystal 
Ball parameters that account for the radiative tail are obtained from the simulation. 

For the ip{2S) — j- J/ip {fi~^ fi~)'7i~^'7i~ decay, both pions are required to have 
Pt > 0.3 GeV/c and the sum of the two-pion transverse momenta is required to be larger 
than 0.8GeV/c. The quantity Q = M( J/^/^Tr+Tr") - M^Ti+n') - M^fi+fi') is required to 
be < 200 MeV/ and to improve the mass resolution the dimuon invariant mass M^+^- 
is constrained in the fit to the nominal J/ip mass value [18]. Finally, both J/ip and ip{2S) 
candidates must have > 2 GeV/ c. The invariant mass spectrum is shown in Fig. [U^b) 
for all selected candidates. For this decay mode the peak is described by the sum of two 
Crystal Ball functions for the signal plus an exponential function for the background. The 
number of signal candidates is 12300±200, the mass resolution is 2.10±0.07MeV/c2, and 
the Crystal Ball tail parameters are fixed to the values obtained from the simulation. 

The fits are repeated in each ip{2S) pt bin to obtain the number of signal and back- 
ground candidates for both decays. 




Figure 1: Invariant mass distribution for all ip{2S) candidates passing the selection cuts 
for the fi~^fi~ decay (a) and the J/ip{^'^ ^~)ti'^tt~ decay (b). 
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4 Cross-section measurement 



The differential cross-section for the inclusive 'ip{2S) meson production is computed from 



da 
dpT 



A^sig(pT) 



jC etot(PT) B ApT 



(1) 



where da/dp^ is the average cross-section in the given bin, integrated over the rapidity 
range 2 < y < 4.5, A''sig(pT) is the number of signal candidates determined from the mass 
fit for the decay under study, etot(PT) is the total detection efficiency including acceptance 
and trigger effects, B denotes the relevant branching fraction and Apx is the bin size. All 
branching fractions are taken from Ref. [IS]: B{ip{2S) e+e") = (7.72 ± 0.17) x 10"^ 
B{ip{2S) J/ipTT+TT-) = (33.6 ±0.4) x lO'^ and B{J/ij /i+/i") = (5.93 ±0.06) x lO'^. 
Assuming lepton universality, we use the dielectron branching fraction B{i(j{2S) — )■ e~^e~) 
in Eq. ([1]), since B{ip{2S) — j- /i'^'yU") is less precisely known. C is the integrated luminosity, 
which is calibrated using both Van der Meer scans [T9|l20] and a beam-profile method |21] . 
A detailed description of the two methods is given in Ref. [22]. The knowledge of the 
absolute luminosity scale is used to calibrate the number of tracks in the vertex detector, 
which is found to be stable throughout the data taking period and can therefore be 
used to monitor the instantaneous luminosity of the entire data sample. The integrated 
luminosity of the data sample used in this analysis is determined to be 35.9 pb~^. 

The total efficiency, etotiPT), is a product of three contributions: the geometrical 
acceptance, the combined detection, reconstruction and selection efficiency, and the trigger 
efficiency. Each contribution has been determined using simulated events for the two decay 
channels. In order to evaluate the trigger efficiency, the trigger selection algorithms used 
during data taking are applied to the simulation. 



55 0.5 

V) 
CM 

% 0.4 
^ 0.3 
0.2 
0.1 



LHCb 
Simulation 



(a) 



2 4 6 8 10 12 14 16 

(GeV/c) 



55 0.2 
to 

^ 0.15 

o 

uT 

0.1 



0.05- 



I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I 



LHCb 
Simulation 



(b) 







2 4 6 8 10 12 14 16 

(GeV/c) 



Figure 2: Total efficiency vs. pt computed from simulation for unpolarized ijj{2S) mesons 
for iP{2S) -> (a) and ij{2S) J/ip{fi+ii-)n+Ti- (b). 



The total efficiency vs. px for the two channels, assuming the ip{2S) meson unpolar- 
ized, is shown in Fig. [21 Extensive studies on dimuon decays of prompt J/ip [H], ip{2S) 
and T [23] mesons have shown that the total efficiency in the LHCb detector depends 
strongly on the initial polarization state of the vector meson. This effect is absent for 
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'ip{2S) mesons coming from 6-hadron decays. In fact for these events the natural polariza- 
tion axis is the ip{2S) meson flight direction in the fe-hadron rest frame, while the ip{2S) 
meson appears unpolarized along its flight direction in the laboratory. Simulations [8] 
and measurements from CDF [2l] confirm this. We do not measure the ip{2S) meson 
polarization but we assign a systematic uncertainty to the unpolarized efficiencies in the 
case of prompt production. Events are generated with polarizations corresponding to the 
two extreme cases of fully transverse or fully longitudinal polarization and the efficiency is 
re-evaluated. The difference between these results and those with the unpolarized sample 
is taken as an estimate of the systematic uncertainty. 

A similar effect exists for the J/ip meson emitted in the ipi^S) J/ip{fi~^fi^)'iT^7r~ 
decay. However, in this case, the ip{2S) meson polarization is fully transferred to the J/ip 
meson since, as measured by the BES collaboration [25], the two pions are predominantly 
in the S'-wave configuratiorQl and the dipion- J/^p system is also in a S-wave configuration. 
This has been verified with data and is correctly reproduced by the simulation. Therefore 
the systematics due to polarization are fully correlated between the two channels and 
we use the systematic uncertainties computed for ip{2S) — )■ fj,~^f^~ also for the il:{2S) — )■ 
J/ipTX^Ti' decay. 

In order to separate prompt ip{2S) mesons from those produced in 6-hadron decays, 
we use the pseudo-decay-time variable defined as t = dz{M/pz), where dz is the separation 
along the beam axis between the ip{2S) decay vertex and the primary vertex, M is the 
nominal mass of the ip{2S) and Pz is the component of its momentum along the beam axis. 
In case of multiple primary vertices reconstructed in the same event, that which minimises 
\dz\ has been chosen. The prompt component is distributed as a Gaussian function around 
t = 0, with width corresponding to the experimental resolution, while for the ip{2S) from 
6-hadron decays the t variable is distributed according to an approximately exponential 
decay law, smeared in the fit with the experimental resolution. The choice of taking the 
primary vertex which minimises \dz \ could in principle introduce a background component 
in the pseudo-decay-time distribution arising from the association of the ip{2S) vertex to 
a wrong primary vertex. The effect of such background is found to be of the order of 
0.5% in the region around t = and has been neglected. The function used to fit the t 
distribution in each pt bin is 

Fit- /p, a, n) = iVsi, 

where N^ig and A'^bkg are respectively the numbers of signal and background candidates 
obtained from the mass fit. The fit parameters are the prompt fraction, /p, the standard 
deviation of the Gaussian resolution function, a, and the lifetime describing the long- 
lived component of ip{2S) mesons coming from 6-hadron decays, r^. In principle, all 
fit parameters are dependent on pt- The function fhkg{t] ©) models the background 
component in the distribution and is defined as the sum of a 5 function and a Gaussian 
function for the prompt background, plus two exponential functions for the positive tail 
and one exponential function for the negative tail, all convolved with a Gaussian function 

^The small fraction of D-wave measured in Ref. [35] has a negligible impact on our conclusion. 
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to account for the detector resolution. The array of parameters is determined from a 
fit to the t distribution of the events in the mass sidebands. 
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Figure 3: Pseudo-decay-time distribution for ^^{23) — /i^yU" (a) and ip{2S) — )■ 
J/ipTT'^'iT^ih) in the pt range 4 < px < 5 GeV/c, showing the background and prompt 
contributions. 



As an example, the pseudo-decay-time distributions for ^'(25') fi^fi~ and ^^{23) 
J/ip TT'^7r~ in the range 4 < < 5 GeV/ c are presented in Fig. [31 The contributions of 
background and prompt ip{2S) mesons are also shown. The values of the prompt fraction, 
/p vs. pt in the rapidity range 2 < y < 4.5, obtained for the and the J/ipn^n~ 

modes, are in good agreement as shown in Fig. HI 
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Figure 4: Fraction of prompt ip(2S) as a function of p^ for the fi~^fi~ mode (solid squares) 
and the J/ipir^iT^ mode (open squares). Error bars include the statistical uncertainties 
and the systematic uncertainties due to the fitting procedure. 



6 



5 Systematic uncertainties on the cross-section mea- 
surement 



A variety of sources of systematic uncertainties affecting the cross-section measurement 
were taken into account and are summarised in Table [TJ 

Table 1: Systematic uncertainties included in the measurement of the cross-section. Un- 
certainties labelled with a are correlated between the fi'^fi~ and J/ipn'^TT" mode, while b 
indicates a correlation between ip(2S) — fi'^fi~ and the J/ip — > fi'^fi" uncertainties [S]. 



Uncertainty source 




J/^ 71+71- {%) 


Luminosity"''' 


3.5 


3.5 


Size of simulation sample 


0.4-2.2 


0.6-1.0 


Trigger efficiency" 


1-8 


1-7 


Global event cuts"''' 


2.1 


2.1 


Muon identification"''' 


1.1 


1.1 


Hadron identification 




0.5 


Track x^"''' 


1 


2 


Tracking efficiency" 


3.5 


7.3 


Vertex fit'' 


0.8 


1.3 


Unknown polarization" 


15-26 


15-26 


Mass fit function 


1.1 


0.5 


Pseudo-decay-time fits 


2.7 


2.7 


-B(V'(25) ^ e+e") 


2.2 




B{^p{2S) J/i) Ti+Ti-) 




1.2 


B{J/ip ^ 




1.0 



A thorough analysis of the luminosity scans yields consistent results for the absolute 
luminosity scale with a precision of 3.5% |22], this value being assigned as a systematic 
uncertainty. The statistical uncertainties from the finite number of simulated events on the 
efficiencies are included as a source of systematic uncertainty; this uncertainty varies from 
0.4 to 2.2% for the ^+ ^~ mode and from 0.6 to 1% for the J/ijj7r+7T~ mode. In addition, 
we assign a systematic uncertainty in order to account for the difference between the 
trigger efficiency evaluated on data by means of an unbiased n+fj^~ sample, and the trigger 
efficiency computed from the simulation. This results in a bin-dependent uncertainty up 
to 8% for the fi+fi~ mode and up to 7% for the J/ip7i+7i~ mode. This uncertainty is 
fully correlated between the two decay modes in the overlapping region. Finally, the 
statistical uncertainty on the global event cuts efficiency (2.1% for both modes) is taken 
as an additional systematic uncertainty [8]. 

To assess possible systematic differences in the acceptance between data and simula- 
tion for the J/ijj 7r+7r~ mode, we have studied the dipion mass distribution. The LHCb 
simulation is based on the Voloshin-Zakharov model [26] which uses a single phenomeno- 
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logical parameter A 



da 



(3) 



dm. 




where $(7??.^^) is a phase space factor (see e.g. Ref. [25]) and in the simulation A = 4 is 

assumed. The dipion mass distribution obtained from the data is shown in Fig. [51 We 
obtain A = 4.46 ± 0.07(stat) ± 0.18(syst), from which we estimate a negligible systematic 
effect on the acceptance (0.25%). Our result is also in good agreement with the BES value 
A = 4.36 ± 0.06(stat) ± 0.17(syst) [25]. 




Figure 5: Dipion mass spectrum for the ipi^S) — t- J/ip-K^ii decay. The curve shows the 
result of the fit with Eq. (|3]) corrected for the acceptance. 

To cross-check and assign a systematic uncertainty to the determination of the muon 
identification efficiency from simulation, the single track muon identification efficiency has 
been measured on data using a tag-and-probe method [27]. This gives a correction factor 
for the dimuon of 1.025±0.011, which we apply to the simulation efficiencies. The 1.1% 
uncertainty on the correction factor is used as systematic uncertainty. The efficiency of the 
selection requirement on the dipion identification has been studied on data and simulation 
and a difference of 1% has been measured between the two. Therefore, the simulation 
efficiencies are corrected for this difference and an additional systematic uncertainty of 
0.5% is included. 

The ip{2S) selection also includes a requirement on the track fit quality. The relative 
difference between the efficiency of this requirement in simulation and data is taken as 
a systematic uncertainty, resulting in an uncertainty of 0.5% per track. Tracking studies 
show that the ratio of the track- finding efficiencies between data and simulation is 1.09 for 
the fi'^fi~ mode and 1.06 for the J/^/^tt+tt" mode, with an uncertainty of 3.5% and 7.3% 
respectively; the simulation efficiencies are corrected accordingly and the corresponding 
systematic uncertainties are included. 

For the requirement on the secondary vertex fit quality, a relative difference of 1.6% 
for the mode and 2.6% for the J/tpir^Ti" mode has been measured between data 





(MeV/c') 
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and simulation. The simulation efficiency is therefore corrected for this difference and a 
corresponding systematic uncertainty of 0.8% {fi'^fi~) and 1.3% {J/ipn^n^) is assigned. 

The systematic uncertainty due to the unknown polarization is computed as discussed 
in Section m The study done for the two extreme polarization hypotheses gives an average 
systematic uncertainty between 15% and 26% for both modes, relative to the hypothesis 
of zero polarization, depending on the pt bin. These errors are fully correlated between 
the two decay modes and strongly asymmetric since the variations of the efficiency are of 
different magnitude for transverse and longitudinal polarizations. 

A systematic uncertainty from the fitting procedure has been estimated from the 
relative difference between the overall number of signal ip{2S) and the number of signal 
candidates obtained by summing the results of the fits in the individual bins. A total 
systematic uncertainty of 1.1% for the fi~^fi~ mode and 0.5% for the J/ipn^TT' mode is 
assigned. 

Finally, to evaluate the systematic uncertainty on the prompt fraction from the ip{2S) 
pseudo-decay-time fit we recompute /p with Tb (see Eq. ([2])) fixed to the largest and 
smallest value obtained in the p^-hin fits. The relative variation is at most 2.7% and this 
value is assigned as a systematic uncertainty on /p. 



6 Cross-section results 

The differential cross- sect ions for prompt ijj{2S) and ip{2S) mesons from 6-hadron decays 
are shown in Fig. [6l where we compare the results obtained for the ip{2S) — )■ and 
ip(2S) — J/ipTi^Ti" channels separately for the prompt and 6-hadron decay components. 
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Figure 6: Comparison of the differential cross-sections measured for prompt ip{2S) (cir- 
cles) and for ip{2S) from 6-hadron decay (squares) in the ip{2S) — > (solid symbols) 
and ip(2S) — > J/ip 7r~^7r~ (open symbols) modes. Only the uncorrelated uncertainties are 
shown. 
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The values for the two cross-sections estimated using the different decay modes are 
consistent within 0.5 cr. A weighted average of the two measurements is performed to 
extract the final result listed in Table |2j 

Table 2: Cross-section values for prompt ip{2S) and ipi^S) from fe-hadrons in different 
Pt bins and in the range 2 < y < 4.5, evaluated as the weighted average of the fi~^^~ 
and J/ip 7r~^7r~ channels. The first error is statistical, the second error is systematic, and 
the last asymmetric uncertainty is due to the unknown polarization of the prompt ip{2S) 
meson. 



Pt [GeV/c] 


rfoprompt r nb 1 


dcTb r nb l 


dpT l-GcV/c-l 


dpT l-GeV/cJ 


0-1 


188 ± 6 ± 18+^2 


22 ± 2 ± 2 


1-2 


387 ± 8 ± 37^^?9 


62 ± 3 ± 6 


2-3 


317 ± 7 ± 26+^^ 


53 ± 2 ± 4 


3-4 


224 ± 6 ± 24+^^ 


39 ± 2 ± 4 


4-5 


135 ± 4 ± 13+^[; 


29 ± 1 ± 3 


5-6 


77 ± 2 ± 7% 


18 ± 1 ± 2 


6-7 


46 ± 1 ± 4+^0 


10 ± 1 ± 1 


7-8 


25 ± 1 ± 2+^ 


6.3 ± 0.4 ± 0.5 


8-9 


14 ± 1 ± 1+2 


3.9 ± 0.3 ± 0.3 


9-10 


8.3 ± 0.4 ± 0.7^?:? 


2.5 ± 0.2 ± 0.2 


10-12 


4.3 ± 0.3 ± 0.4+°| 


1.4 ± 0.1 ± 0.1 


12-16 


1.5 ± 0.1 ± 0.2+°;^ 


0.51 ± 0.04 ± 0.06 



The differential cross-section for promptly produced ip{2S) mesons, along with a com- 
parison with some recent theory predictions [281431] tuned to the LHCb acceptance, is 
shown in Fig. [71 In Ref. |2B] and Ref. [21] the differential prompt cross-section has been 
computed up to NLO terms in nonrelativistic QCD, including colour-singlet and colour- 
octet contributions. In Ref. [301 [31] the prompt cross-section has been evaluated in a 
colour-singlet framework, including up to the dominant NNLO terms. Experimen- 
tally the large-pT tail behaves like p^^ with /3 = 4.2 ± 0.6 and is rather well reproduced, 
especially in the colour-octet models. 

The differential cross-section for ^/'(2S') produced in 6-hadron decays and the compari- 
son with a recent theory prediction [32] based on the FONLL approach [6l[7] are presented 
in Fig. [HI The theoretical prediction of Ref. [32] uses as input the h — )■ il){2S)X branching 
fraction obtained in the following section. Experimentally the ^'(25') mesons resulting 
from 6-hadron decay have a slightly harder spectrum than those produced promptly: 
/3 = 3.6 ± 0.5. By integrating the differential cross-section for prompt ip{2S) and 'iIj{2S) 
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from 6-hadrons in the range 2 < y < 4.5 and <16 GeV/c, we obtain 

fTprompt(^(2^)) = 1.44 ±0.01 (stat) ±0.12 (syst)+°:20 (pol) ^ib, 
ab{'4^{2S)) = 0.25 ± 0.01 (stat) ± 0.02 (syst) yih, 

where the systematic uncertainty includes all the sources listed in Table [H except for the 
polarization, while the last asymmetric uncertainty is due to the effect of the unknown 
ip{2S) polarization and applies only to the prompt component. 
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Figure 7: Differential production cross-section vs. for prompt ^(25*). The predictions 
of three nonrelativistic QCD models are also shown for comparison. MWC [28] and 
KB [29] are NLO calculations including colour-singlet and colour-octet contributions. 
AL [30|[3T] is a colour-singlet model including the dominant NNLO terms. 



7 Inclusive b ip{2S)X branching fraction mea- 
surement 

The inclusive branching fraction for a 6-hadron decaying to ijj{2S) is presently known with 
50% precision: B{b — )■ ip{2S)X) = (4.8 ± 2.4) xlO~^ [18]. Combining the present result for 
ab{'ip{2S)) with the previous measurement of ab{J/ip) [H] we can obtain an improved value 
of the aforementioned branching fraction. To achieve this, it is necessary to extrapolate 
the two measurements to the full phase space. The extrapolation factors for the two 
decays have been determined using the LHCb simulation [12] and they have been found 
to be a4,r(<//V')=5.88 [8] and a4,r('^(25'))=5.48. Most of the theoretical uncertainties are 
expected to cancel in the ratio of the two factors ^ = a4-„{ip{2S)) /a^-^^J/ip) = 0.932, 
which is used in Eq. (jlj. A systematic uncertainty of 3.4% is estimated for this correction 
and included in the final result below. Therefore 

Bjb ^ ^i2S)X) ^ abmS)) 
B{b^J/^X) ^ ab{J/ij) ■ ^ ' 
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Figure 8: Differential production cross-section vs. for ip{2S) from 6-hadrons. Tfie 
shaded band is the prediction of a FONLL calculation [S1EIE2]- 



For (Th{J/ip) we rescale the value in [8] for the new determination of the integrated 
luminosity {C = 5.49 ± 0.19 pb~^). For ab{ip{2S)) we use only the data from the 
ip{2S) — fi'^fi~ mode to cancel most of the systematic uncertainties in the ratio. Effects 
due to polarization are negligible for mesons resulting from 6-hadron decay. We obtain 

B(b ^ iIj(2S)X) , , , , 

— ' ' = 0.235 ± 0.005 (stat) ± 0.015 (syst), 

where the correlated uncertainties (Table [1]) between the two cross-sections are excluded. 
By inserting the value B{h J/i^X) = (1.16 ± 0.10) x 10"^ [TB] we get 

B{b i^{2S)X) = (2.73 ± 0.06 (stat) ± 0.16 (syst) ± 0.24 (BF)) x 10~^ 

where the last uncertainty originates from the uncertainty of the branching fractions 
B{b J/iIjX), B{'iIj{2S) -> e+e") and B{J/ij 

The ratio of the ip{2S) jX^ to J/ip — jX^ jjT differential cross-sections is shown vs. 
Pt in Fig. [9] for prompt production (i?p, Fig. M^a)) and when the vector mesons originate 
from 6-hadron decays {Rb, Fig. [9]^b)). Since it is not known if the promptly produced 
'0(25') and J/ip have similar polarizations [33], we do not assume any correlation of the 
polarization uncertainties when computing the uncertainties on Rp. The increase of Rp(b) 
with is similar to that measured in the central rapidity region by the CDF [23] and 
CMS [i] collaborations. 

8 Conclusions 

We have measured the differential cross-section for the process pp — )■ tlj{2S)X at the 
centre-of-mass energy of 7 TeV, as a function of the transverse momentum in the range 
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Figure 9: Ratio of ^^{23) — > fi~^fi^ to J/ip — )■ fi~^fi^ cross-sections for prompt 
production (a) and for 6-hadron decay (b), as a function of p^. 



Pt{^{2S)) < 16 GeV/c and 2 < y{^{2S)) < 4.5, via the decay modes ^^(25) ^ and 
ip{2S) — > J/ifjix^n" . The data sample corresponds to about 36 pb~^ collected by the LHCb 
experiment at the LHC. Results from the two decay modes agree. The ip{2S) prompt 
cross-section has been separated from the cross-section of ip{2S) from 6-hadrons through 
the study of the pseudo-decay-time and the two measurements have been averaged. In 
the above kinematic range we measure 

cTprompt(^(25)) = 1.44 ± 0.01 (stat) ± 0.12 (syst)l°;^° (pol) ^ib, 
ab{il){2S)) = 0.25 ± 0.01 (stat) ± 0.02 (syst) ^ib. 

The measured ip{2S) production cross-sections are in good agreement with the results of 
several recent NRQCD calculations. In addition, we obtain an improved value for the 
b — tlj{2S)X branching fraction by combining the two LHCb production cross-section 
measurements of the two vector mesons J/ip and ip{2S) from 6-hadrons. The result, 

B{b 'ipi2S)X) = (2.73 ± 0.06 (stat) ± 0.16 (syst) ± 0.24 (BF)) x 10-^ 

is in good agreement with recent results from the CMS collaboration [9] and is a significant 
improvement over the present PDG average [TSj . 
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